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The exact Bethe-Salpeter amplitude for the fermion-antifermion bound state in the Schwinger 
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I. INTRODUCTION 

The formation of bound states is one of the most fun- 
damental unsolved problems in quantum field theory. It 
becomes especially troublesome if relativistic nature of 
these states cannot be neglected and retardation effects 
in the interactions among constituent particles have to 
be taken into account. 

In the fifties of the last century, the quantum field the- 
oretical equation for the bound state amplitude - the 
so called Bethe-Salpeter (B-S) equation - was formu- 
lated [2 [2] and the framework for considering bound 
states was established. However, for serious reasons, this 
equation cannot be solved in realistic quantum field the- 
ory. From the mathematical point of view, it is an ex- 
tremely complicated multidimensional integral equation 
for a kind of a wave-function but the main difficulty 
conies from the fact, that from the very beginning one 
needs to know the nonperturbative propagators for in- 
gredient particles and their interaction kernel. Obviously 
none of these quantities is known in theories like QED or 
QCD or even in simpler field theoretical models. 

The significant feature of the relativistic B-S amplitude 
is its dependence on the so called relative time, or rela- 
tive energy, the role of which in characterizing a bound 
state has not yet become completely clear. Their pres- 
ence in the theory mirrors the effects of retardation in 
the interactions and constitutes the main difference with 
respect to the nonrelativistic studies. It is also a source 
of mathematical (and interpretational) complications. 

For the above reasons over the last sixty years much 
work has been devoted to the investigation of approx- 
imated versions of the equation [3l H]. These approx- 
imations consisted on simplifying the propagators and 
interaction kernel (the most often the 'ladder' approxi- 
mation was used) and of reducing the equation to lower 
dimensions by neglecting the relative time or relative en- 
ergy dependence [5HTU]. In this case retardation effects 
and consequently a part of physical information were lost 
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not allowing for the full study of the bound state prop- 
erties as for instance its threshold structure. Some in- 
sight into the relative time (or energy) dependence of the 
bound system has been obtained in certain models |111 - 
I14j . The wave-function dependence of relative variables 
was mainly given in euclidean space |15H21| . where it is 
easier to obtain, but the Minkowski-space approach for a 
bound system has also been developed pTH^ . 

A very unpleasant situation is that full results are still 
lacking even in simple field theoretical models. There 
are only few examples, as for instance the so-called Wick- 
Cutkosky Model [1211301131]: where it was possible to find 
'exact' solutions (but with the 'ladder' approximation 
embodied into the model). In this context we would like 
to mention the result obtained in the Schwinger Model 
(SM) - a two-dimensional massless electrodynamics |31] 
- where full (i.e. exact) B-S wave-function was found, in 
the spirit of [33l [S^ , from the residue of the pole corre- 
sponding to the bound particle, but also as a direct com- 
putation of the field amplitude [351 I3S] opening thereby 
the possibility of its investigation. 

To our knowledge, no other truly exact solutions of 
the B-S equation in nontrivial field theoretical models 
are known. Our main motivation in this work is then 
the examination of the previously found B-S amplitude 
from the point of view of the relative time and energy 
dependence. One should stress here that the considered 
model is fully relativistic, since it deals with massless 
fermions. 

Although B-S function docs not have any direct prob- 
abilistic interpretation, it appears in the matrix elements 
and scattering processes involving bound states [311 137] ■ 
The results may be then of certain interest both in parti- 
cle and nuclear physics. As pointed out in [llj the relative 
time or energy dependence is important for the processes 
in which the bound state appears in loop integrations, 
but also can manifest itself when one deals with high 
momentum transfers. The results, which in the SM are 
exact, may constitute also the measure of accuracy of 
various approximation schemes. 

The Schwinger Model offers one more attractive possi- 
bility. As it is well known, it is the theory with nontrivial 
vacuum [3HH131 , which has the character of the 0- vacuum. 



Apart from the confinement of fermions this feature re- 
sulting from the the existence of instanton sectors makes 
the theory to be a simple model of QCD. The nonzero in- 
stanton sectors do contribute to the B-S amplitude, and 
the appropriate contribution was explicitly found in |35j . 
One can then analyze how instantons influence the rel- 
ative energy and momentum dependence of the bound 
state wave-function (an attempt to account for instanton 
effects in bound states was made in [HI 03]). 

The present paper is organized as follows. In the next 
section we invoke the exact B-S amplitude in coordinate 
space for all contributing instanton sectors. We inves- 
tigate its position dependence for fixed values of rela- 
tive time and vice versa in the center-of-mass frame and 
present the appropriate plots. In section |III| we concen- 
trate on the behavior of the amplitude in the relative 
momentum for certain characteristic values of relative 
energy again in the center-of-mass frame and we depict 
the results of the numerical computations for each sec- 
tor. At the end we present the density plots revealing 
the dependence of the wave-function on the whole two- 
momentum argument and show the trajectories of the 
threshold cusps. And finally in the last section we give 
the summary of the obtained results and propose some 
conclusions. 



II. BETHE-SALPETER AMPLITUDE IN 
COORDINATE SPACE 

In order to fix the notation we start this section with 
recalling basic facts concerning SM. Its Lagrangian den- 
sity is given by 

C{x) - *(.x) (irO^ - gA'^ixh^^'i'ix) 

-V''(x)F^.(x)-^(a^^^(a:)f, (1) 

where g is the coupling constant and A the gauge fix- 
ing parameter (in what follows we choose Landau gauge 
setting A — >■ c»). Dirac gamma matrices are chosen in 
two-dimensional form: 
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7^ = 7°7^ 
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(2) 



and for the metric tensor we have 0°° = —0^^ = 1. In the 
Schwinger's papers |32j the nonperturbative formula for 
the fermion propagator in coordinate space was found. 
It may be given the following form: 



S{x) = So{x)exp [~ig'^/3{x)] , 
with Sq{x) being the free propagator: 

1 ^ 



So{x) = - 



27r x^ — is 



(3) 



(4) 



and function /3 defined by 

/?(x) = (5) 

-f + 7B + In v/;?^V4 + fH^o'\y/J^^)] : 



2g2 



2ff^ 



timelike, 



7£ + In v/=Ai2^74 + Koi^/^i^:^) 
X spacelike. 



Symbol je denotes here the Euler constant, /i is the 
Schwinger boson mass (p"^ = g'^ /it) and H^ ' and Kq 
are Hankel function of the first kind and Basset function 
respectively. 

In the further development of the model higher Green's 
functions in explicit form were obtained. Special atten- 
tion was paid to two- fermion function [46H48| . The anal- 
ysis of its analytical structure showed that it contains the 
bound state pole with the residue defining the B-S am- 
plitude [35l [36] . It has the contributions from fc = and 
fc = ±1 instanton sectors and may be given the form: 



where 



$ 



<^p{x) = '^f{x) + '^^p\x), 



f {x) = -2^/^S[x)"t^ sm{Px/2), 



.(1) 



^'p'{x)-^ 



_H_^lE^^g''p{^)^-^er'^ cos(Px/2) 



(6) 

(7) 
(8) 



P denotes here the total two-momentum of the bound 
state (satisfying P^ = /i^), x = [t,r] is the relative two- 
position variable and 9 is the parameter defining the 9- 
vacuum. 

The bound state amplitude is then defined entirely by 
the function j3 and trigonometric functions. The presence 
of these functions (sine for A: = 0, and cosine for fc = ±1) 
shows that there are regions where noninstantonic sector 
dominates, and other, where instanton contributions are 
larger. 

For our further analysis we choose the center-of-mass 
frame in which P = [/^, 0]. In such a case the trigono- 
metric functions depend only on the relative time and 
expressions nh and ([S]) reveal quasi-periodic structure in 
this variable. For t approaching (2n + l)7r//x the am- 
plitude <i>p dominates over topological part, but for t 
close to 2mT I jj, the instanton corrections become more 
important (except the light cone). 

If we fix the relative time to take certain chosen value, 
the dependence of the amplitude on the spatial variable 
r may be easily plotted. Since $p and $p are 2x2 
complex matrices, we will draw the behavior of the quan- 
tities: 
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FIG. 1: The dependence of amplitudes |"l>p | (dotted line), |$p | (dashed line) and |"I>p| (solid line) in units of fi on relative 
position r with fixed relative time t. The variable t for the successive plots is chosen to be: 0, 1, 3, 6, 9, 12. All coordinates 
are in units of /i~^. The dotted hue for i = 3 and f = 9 is covered by the sohd line since the instanton contribution is tiny. 



without any claim to their probabilistic interpretation. ments: 
With this definition the dependence of |$p | on the 
unphysical ^-parameter (after all in the SM it may be 
gauged away) disappears. The real and imaginary parts 
of (It]) and (Is]) may also be easily plotted if needed. 
Thanks to the trace properties of Dirac gamma matri- 
ces, for the total amplitude, we get 



l$p| = 



i tr[$+$p] = 



4°)P 



|$wp. 



(11) 



In Figure |T]th c dependence of both contributions to- 
gether with (11) on the relative position r for fixed rel- 
ative time is plotted. For i — > the amplitude \^p \ 
(which is indicated in all graphs with a dotted line) is 
hardly visible, since it becomes a Dirac delta function, 
which is not graceful for drawing. Due to the denomina- 
tor in ^ the amplitude for fc = exhibits the singulari- 
ties on the light cone. This does not refer to the nonzero 
instanton sector (dashed line) . The solid line corresponds 
to the expression \^p\. 

On the third plot (i.e. for t = 3 in units of fi^^) we ob- 
serve the broadening of peaks in \^p | and strong reduc- 
tion of the value of |^p |- This is understandable since 
3 is close to tt (hence sin(/ii/2) « 1 and cos(/it/2) « 0). 
The same refers to t = 9, which approaches Stt. On the 
other hand, for i = 6 (i.e. close to 2tt) and t = 12 (i.e. 



s(i) 



(0), 



almost 47r) |$p | becomes more important, and |$p | is 
close to zero, apart from the proximity of the light cone, 
where peaks turn now very narrow. This kind of light- 
cone singularities is connected with the masslessness of 
the fundamental fermions in the model. 

Since the Basset function tends to zero for large argu- 



Koiz) 



it may be easily seen from (l7| and (|8| that both ampli- 
tudes behave like 1/ Vrf ^^'^ become comparable to each 
other for large separations. This is a relatively slow de- 
crease and is again related to the formation of the massive 
bound state by massless fermions and to the presence of 
fermionic zero modes in the instanton background. One 
should also note, that in that region (strictly speaking 
outside the light cone) the whole amplitude in both sec- 
tors is real (for 9 = 0). 

If we now fix the relative position and investigate in an 
analogous way the relative time dependence of the am- 
plitudes, we expect their quasi-periodic behavior due to 
to the presence of trigononretric functions (in a boosted 
frame, this character would appear also in the r depen- 
dence). They make \^p' \ to have (by turn) nodes at 
the integer values of tt/ fi. This can be seen in Figure p] 
together with the typical light-cone singularities. Again 
the peaks of |$p | widen for odd values of n/fi (what is 
also in general accompanied by the damping of instan- 
ton contributions) and become sharp for even ones. The 
latter are regions of instanton domination. On the other 
hand |$p| is much more smooth with no nodes at all. 
For example if the relative distance gets large, there is 
for relatively small t a whole range of almost constant 
|$p|. One can easily check, that in this case: 



|4°)(a;)|: 



!LplB/2^Ko{y/^J^^)\ 



sin(Fa;/2)|, (12) 



\<^^p\x)\ = J^e^^/^e-^'^°'-^^^^^\cosiPx/2)\. 
V 27r 




r=6 




r=3 



-3 n —7t 




n 3 TT 



r=9 



-3n -IT n "in 




- t 




-3 n —n 



K 3 JI 



-3 n -n n 3 n 



FIG. 2: Same as Figure [T] but with fixed relative position r. The variable r for the successive plots is chosen to be: 0, 1, 3, 6, 
9, 12. 



Since for those fairly large values of 



the function 



K^^lyJ—jjflx^) equals almost zero, then both contribu- 
tions complement each other and |$p| becomes a very 
slowly varying function. This interplay between oscillat- 
ing noninstantonic and instantonic parts constitutes an 
interesting observation. 

For large values of relative time (i.e. for large timelike 
separations) we have to use the upper formula of ([5]). The 
asymptotic behavior is then dictated by the presence of 
a logarithm since for the Hankel function we have: 



H'^\^) 







and we conclude that the amplitudes fade out as l/v^, 
hence similarly to the previous case. 



III. BETHE-SALPETER AMPLITUDE IN 
MOMENTUM SPACE 

A. fc = instanton sector 

In real four-dimensional quantum field theory calcu- 
lations one mostly has to do with Feynman diagrams in 
momentum space. In this representation the 5-matrix el- 
ements have their natural form and the analytical prop- 
erties of the Green's functions (poles, cuts) are closely 
related to the physical quantities. One can expect that 
also B-S wave-function should reveal certain important 
properties regarding for instance the internal threshold 
structure. 

To analyze the dependence of B-S wave-function on 
the relative two- momentum variable Q = [E^ g], what we 



need, is the transformed propagator ([3]), since 

<^f{Q)^-2V^ j cPxe'^''sin{Px/2)S{x)-i'^ (13) 
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Due to its complicated form it cannot be, however, ex- 
plicitly Fourier-transformed. The direct numerical trans- 
formation is, on the other hand, very slowly convergent 
due to the oscillatory nature of integrand function. For 
these technical reasons, we prefer to make use of the inte- 
gral representation of the euclidean version of propagator 
S{p) |49| . which may be in the straightforward way con- 
tinued to spacelike momenta: 



S{p)^ 



^ 



p7e/2 



2 (-p2)5/4 

dx\/xJi{x) exp 



(14) 



where Ji is the Bessel function. One might say that 
in this formula the angular part of the two-dimensional 
transform has been performed, and the radial part is still 
left. To make the appropriate plots we will need to know 
the values of S{p) in the who le momentum plane. The 



analytical continuation of (14) to the region, where p > 
0, gives: 



S{p)-- 



(l_i)eT«/2 



2 (p2)5/4 

: / dx-^a; Ji (x) exp — -ffg (m^/vP^) 



(15) 



Using together (14 1 and (15) we get the integral repre- 
sentation of (13). 



Figure p^ shows the dependence of the quantity |<i>p | 
defined again by the formula ([9]), (but referring now to 
the object ( jTs] )) on the relative momentum for two chosen 
values of relative energy: E = and E ~ 1/A ^ \n the 
center-of-mass frame. 
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FIG. 3: The dependence of the amplitude |$p | (in units of 
fi^^) on the relative momentum q for fixed values of relative 
energy {in ^) : E — (upper plot) and E = 1/4 (lower plot). 



The observed singularities are of kineniatical nature 
(similar ones appear in [5BH25]) and correspond to the 
external 'legs' of the B-S wave-function ($p is not an 
amputated function). They arise, when (P/2 + Q)^ — 
or {—P/2+Q)'^ = 0. These two conditions may be rewrit- 
ten as 



{n/2 + E)^ - q^ = 0, 
{^^i/2 + Ef-q^ = 0. 



(16) 



Hence for _E = the singularities are situated at g = 
±1/2 /x. If relative energy increases, four such singu- 
larities appear. For instance a.t E = 1/4 /i, we get 
q = ±1/4 n and q = ±3/4 /i. This is clearly visible 
on the second plot of Figure |3] 

When the relative energy reaches the threshold corre- 
sponding to the Schwinger boson, a new effect appears. 
In general n-th order thresholds result from satisfying the 
equations: 



(P/2 + Qf = {npif 
{-P/2 + Qf = (nM)^ 



(17) 



i.e. when appropriately high momentum transfer be- 
comes possible. In the center-of-mass frame their appear- 
ance requires that (±/i/2 + EY — n^/i^ > and hence: 



E^ 



nT 



1 



E 



n± 



1 



>0. 



(18) 



For n = 1 we see that, while taking upper signs in (18 1, 
the first two cusps should appear, \i E > 1/2 /i. For 
instance as E = 3/4 /^, the solutions of the first equation 



of (17) are 




(19) 



These cusps are clearly visible on the first plot of Fig- 
ure [4] in places exactly predicted, apart from kinematic 
singularities at g = ±1/4 /i and q = ±5/4 /^ resulting 
from (fiel. 
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FIG. 4: Same as in Figure [s] but for E = 3/4 (upper plot) 
and E = 7/4 (lower plot). 

Moving further we see, that if relative energy exceeds 
the value of 3/2 /i, the condition ( [l8| with lower signs 
may be satisfied too. This should produce additional 
two cusps. For exemplary value E — 7/4 fi it may be 
easily found that the positions of all four cusps are 



,3 



65 
g==±^-— ^ w±2.02 ij,. 



This is exactly, what we see on the second plot of Fig- 
ure E^ (beside kinematic singularities which now moved to 
±5/4 /x and ±9/4 fi). 

Putting n = 2 into (lls]) we find that second order 
cusps connected with two-boson thresholds appear for 
E > 3/2 fi (the first two) and for E > 5/2 fi (the other 
two). These cusps are however extremely weak (they are 
smaller than the thickness of the line we use), although 
visible on very high resolution plots. For instance on 
the lower plot of Figure |4J the n = 2 cusps appear at 
q — ±vT7/4 fi sa ±1.03 fi. We have to stress that we pos- 
itively verified also the appearance of all four new cusps 
for E = 11/4 fi (they are at ±7105/4 fi and ±\/T7/4 fi) 
but they are so faint, that presenting the appropriate 
drawings here is useless. 



B. 



fc = ±1 instanton sectors 



To find out whether instanton terms exhibit the same 
threshold structure we have now to Fourier transform the 
expression ([8]), which, however, similarly as for fc = 0, 
cannot be done explicitly. Fortunately an analogous in- 
tegral representation for Q^ < as in (14) may be used. 
If we define the quantity 



B{p) = J' 



Pxe^'^^e*^ ^(2:) 



(20) 
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FIG. 5: Same as in Figure [s] but for j$<,^' 



we obtain 



B{p)^-2m 



fi (-p2)3/4 



-7e/2 



(21) 



X / (ix-v/a;Jo(a;) exp 
Jo 

where Jq is the Bessel function. 

The instanton contribution to B-S function may be 
now written as 



4V7r 



(22) 



5(f + g)+s(-f 



For the timclikc Q, one finds 



B{p) = - 



2(1 -i) 1 



= -7e/2 



(23) 
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FIG. 6: Same as in Figure [i] but for 1$^^']. 

As may be seen in Figures [5^ and ^ the amplitude |$p | 

has analogous kinematic singularities as |$p |, whose lo- 
cations are dictated by ( [l6| . At the same threshold val- 
ues of E as before (i.e. 1/2 fi, 3/2 fi ior n — 1 and 3/2 fi, 
5/2 fL for n = 2 and so on) and for the same relative mo- 
menta the resonant cusps occur. Their positions are at 
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FIG. 7: The density graphics showing the behavior of |"I>p | (left diagram) and \^p \ (right diagram). Bright areas correspond 
to regions of large values and dark areas of small ones. The colors are relative within each plot so they should not be compared 
between left and right diagram. 



±3/4 /x on the upper plot of Figure|6](£' = 3/4 n) and at 
±3/4 ^, ±765/4 n on the lower one {E = 7/4 n). They 
are now much less pronounced although still perfectly 
visible. 

It is important to see, how these peaks move with rel- 
ative two-momentum. If they really correspond to the 
Schwinger boson thresholds, they should be located on 
hyperbolas ( [l7| in the plane (E,q). That it is actually 
the case may be seen from the pictorial density diagrams 
of Figure [7] where light areas represent large values of the 
amplitudes, and dark small ones. Apart from the white 
lines, which show the locations of kinematic singularities, 
the two hyperbolas corresponding to rt = 1 thresholds 
are noticeable. We checked that the n — 2 hyperbolas 
are also present, but their observation require very pre- 
cise color plots and simultaneous artificial lowering the 
large values of |$p | in order to reduce the maximal 
amplitude of colors. 

The important question arises, how these effects 



emerge in (131 and (22) from the mathematical point 



of view. Why different expressions for |<&p | and |<&p | 
display almost identical threshold structure? The crucial 
observation is that infinite integrals (15) and (23) owe 



their convergence to the presence of oscillating (phase) 
factors coming from Bessel and Hankel functions. If we 
expanded the exponential into a power series, the n-th 
term under the integral would contain a factor (I stands 
for or 1): 



Jiix) 



H, 



(1) 



A* 



(both for / = and ^ = 1) contributes e^*^ to the phase 

factor, and n Hankel functions e'"'^/vP^ ^. If both factors 
cancel, i.e. if 



M 



= ±1, 



(24) 



(of course only '+' comes into play) the value of integral 
over X rapidly increases. The condition ( 24 ) corresponds 
exactly to (17). This again proves the correctness of our 



Now, in the large x asymptotics, the Bessel function 



former observation, that the cusps are localized on hy- 
perbolas in the relative two-momentum plane. 



IV. SUMMARY 

Summarizing the obtained results we would like to 
stress in the first place, that Schwinger Model turned 
out again to be exceptionally useful theory. The possi- 
bility of investigation exact bound-state wave-function in 
relativistic quantum field theory is unique even among 
model studies. Special mention deserves the fact that it 
possesses two features typical to QCD: confinement and 
instantonic vacuum. 

The full Bethe-Salpeter wave-function in coordinate 
space in this model is known from the residue in the 
bound-state pole of the two-fermion Green's function, 
without the necessity of solving B-S equation itself (to 
tell the truth even the interaction kernel needed to for- 
mulate this equation is unknown). The contributions to 
the wave-function come from k — and k — ±1 instan- 
ton sectors (but not higher, since this quantity is defined 
through the operator bilinear in fermion fields). In sec- 
tion In] it was analyzed in the center-of-mass frame, how 



this wave- function behaves in relative position and in rel- 
ative time. It turns out that it displays the light cone 
singularities and decreases as \l \fr for large space-like 
separations and l/vi for time-like ones. This result may 
be, however, typical for D — 2. As functions of relative 
time both noninstantonic and instantonic contributions 
show the oscillatory behavior shifted by a phase of 7r/2 
with respect to each other. Owing to that the total func- 
tion defined by (111 is smoothed and has no nodes. 

In the Fourier space we analyzed the behavior of I'I'p | 

and I'I'p I on relative momentum with fixed relative en- 
ergy. In both instantonic sectors we found kinematic sin- 
gularities coming from the external 'legs' of the unam- 
putated amplitude. For E exceeding the threshold val- 
ues there appear peaks corresponding to the Schwinger 



boson resonances. They emerge both for fc = and 
fc = ±1. They are located on hyperbolas defined by equa- 
tions p7|). This threshold behavior has been suggested 



in approximated scalar-scalar model jllj . 

Since the Schwinger Model is on the one hand a non- 
trivial field theory and on the other is analytically solv- 
able, the obtained results may be of certain importance 
for studies of more advanced QFT models. Both the 
retardation effects and going beyond 'ladder' approxima- 
tion are significant from the point of view of the correct 
description of the bound state [SD] . The model opens also 
the possibility of investigations of the B-S amplitude in 
the boosted frame, the analytical properties in the com- 
plex relative energy plane and of form-factors. 
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